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Abstract 

Basing on the picture of Frenkel excitons in molecular crystals described by 
the XY-model we consider the linear response of linear cyclic aggregates at finite 
temperature. The exact results for characteristics of the response are obtained. 
In particular, we calculate time-dependent two-point correlation functions at finite 
temperature for the cyclic 1-D XY-model. 
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1 Introduction 



Last five years a number of papers were devoted to the consideration of the so-called J- 
aggregates, molecular linear aggregates with redshifted unusually sharp absorption bands 
(for the recent review see Q). Such aggregates give the infrequent possibility to describe a 
nontrivial physical situation by some sort of 1-dimensional quantum Hamiltonians which 
could be treated exactly. To be more precise, the optical excitations of the system (exci- 
tons) after the second quantization procedure could be described by the Hamiltonian of 
spin chains (we will see in section 2 that this is the anysotropic XY-chain). Moreover, 
theoretical investigations of J-aggregates essentially use the information about the struc- 
ture of the exact eigenvalues and eigenfunctions of the Hamiltonian which are well-known 
due to the paper by Lieb, Mattis and Schults since 1961. 

The purpose of this paper is to improve the results on the linear response of the 
J-aggregates in the following aspects: 

1. As we will remind in the next section, there is no freedom to restrict the con- 
sideration from the XY-model to the XX-model, because the relative anysotropy 
parameter is the order of unit. This fact causes a number of complications which 
are not only technical ones. Indeed, there is no more good quantum number like 
the number of particles in the chain (what is obvious from physical features of the 
problem) and one has to deal with a full Fock space of particles in question. 

2. The same problem (i.e. the calculations in a full Fock space of the theory) appears 
even in the XX-model when the effects of nonzero temperature are considered and 
the correlation functions are some traces on a full Fock space. 

Both of these problems are treated in the paper. 

The next question of interest is a choice of the boundary conditions for finite linear 
aggregates. It is well understood that the choice of the general closed boundary conditions 
leads to the new nontrivial problem. Indeed, the exchange statistics for the excitons (which 
are statistically Paulions, i.e. fermions on site and bosons on different sites) is irrelevant 
for the open boundary conditions because no real permutation is allowed. The situation 
changes if the linear aggregate is a circle. It leads to a nontrivial loop which gives the 
possibility to permute particles via the 'glued' boundary. Then the effects of the exchange 
statistics exhibit themselves ||. These effects are negligible in the continuous limit and 
were not investigated when the main interest was concentrated on the thermodynamics 
of the chains. A large number of the classical and recent papers concern the subjects 
B H 0' HI- But ^ ne effect of the closure is relevant and important for finite systems 
like J-aggregates. Moreover, they allow to trace the crossover from small to macroscopic 
chains jl]. Namely the cyclic finite aggregates will be considered in this paper. 

We have to note that the exact treatment of excitons in the aggregates on the base of 
the cyclic XX-model were considered in [[|. In Section 3.4 we will criticize the results of 
Ref.i. 

The paper is organized as follows. In the next section we give a brief description of 
excitons in molecular chains which are basic for understanding of the subject. While 
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opticians began to deal with excitons long time ago, mathematical physicists (who are 
mainly concentrated on exact solutions) are not so familiar with the field. The section 
can be considered as an introduction in the excitonic physics. It also contains the formula 
for the linear response function in terms of correlators of the system in the convenient for 
the following form. Section 3 contains a main new result of the paper. It is devoted to the 
exact calculation of two-point correlation functions for the cyclic XY-model in a closed 
form. The conclusion completes the paper with a discussion of the resulting formulae for 
characteristics of the linear response for the cyclic J-aggregates. Appendix A contains the 
evaluation of the main integrals of Section 3. 



2 Linear response of excitonic system 

In the present work we will calculate the optical response of a gas of molecular aggregates. 
We are interested in the average dipole momentum induced by an external time dependent 
electromagnetic field in the linear approximation with respect to the field amplitude . 

The most general linear relation between the dipole momentum and the field has the 
form 



Pa(t) = J a a p{T)Ep{t - r)dr 
o 



where p a is the dipole momentum per unit volume, Ep denotes the external field, a a p{r) 
is the linear response function and a, (3 are tensor indices. 

We suppose that the aggregates are independent, i.e. we will neglect the interaction 
of the aggregates. In this case the dipole momentum per unit volume is just a product 
of the aggregate concentration v and the dipole momentum of a single aggregate. Hence 
the response a a p is also proportional to the concentration v. 

where cr a p(r) is the response function of a single aggregate. 

As the model for a molecular aggregate we consider the 1-D closed M-sites molecular 
chain (that can be regarded for example as a model for the benzene vapor under a low 
pressure). Let us derive the excitonic Hamiltonian of the aggregate following to the 



original paper | ID| . We start with the Hamiltonian of the molecular aggregate in the 
form: 

M i 

H =E#n+2£V nm (1) 

n=l n,m 

Here H n is the Hamiltonian of a molecule on the n-th site, V nm is the operator the of 
interaction between n-th and m-th sites, J2' means the absence of the term with n = m. 

As far as the molecules are supposed to be electrically neutral, the first approximation 
for the interaction V nm is dipolar one. 

Let us denote the complete set of eigenfunctions of the free molecular Hamiltonian H n 
by {0{} and the corresponding occupation number by N n f . After the second quantisation 
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procedure the state of the aggregate becomes a function of occupation numbers 

|...iV»/...> 

and Hamiltonian (|]) can be rewritten by the following way: 

H = E e / & nA/ + ^ £' bl f ,bl g ,b mg b nf (fg'\V nm \fg) (2) 

n,f n,m 

f,f',g,g' 

The creation (annihilation) operators b^t {b n f) of the /-state on the n-th site satisfies 
the standard commutation relations: 

{bnfMj} = 1 



bnfbnf = b+ f b+ f = 
with the constraint (each molecule can be in a single eigenstate only), 

/ 

e/ is the energy of the f-th excited state. 

Further we shall assume that the lowest excited molecular state plays the main role 
in the absorption and radiation processes and therefore the index / runs only over two 
values (ground state) and / (the first excited state). This allows us to present the 
Hamiltonian in the form: 

H = H 1 + H 2 + H 3 + H 4 (3) 

where the first term describes the energy of nonexcited molecules, 

H 1 = Me o + ~£'(00|Km|00), (4) 

^ n,m 

the second term is the energy of excitation of a single molecule in the aggregate: 

H 2 = E(Ae / +P / )iV n/ 

n 

Ae / = e f ~ e o 

Vf = T! {(o/ \v nm \ of) - (00 \v nm \ oo)} 

m 

and the rest terms have the form: 

n,m 

H 4 = |E ' M lm 1 bnobmobmfbnf + b+fb^fbrnobnO \ 
n,m 1 
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We have neglected here the term 

\t! ((oo Km\ oo) - (0/ \v nm \ of) + (// \v nm \ //)) b+ f b+ f b mf b nf 



n,m 



which arises due to the multipolar interaction and is small comparing with H 3 and H A . 
The matrix elements 

M/ m = (0/|Km]/0>, 

(00 \V nm \ ff) ~ M„' m , (5) 

(ff\V nm \00)~Mf m , (6) 

characterize the excitation 'hopping' from the site n to m. Introducing new (excitonic) 
operators 

Pn = Kobnf , Pn = Kf^nO 

which play the role of the creation (P„) an d annihilation (P n ) operators of the exciton 
on the n-th site, one can write the Hamiltonian in the finite form: 

H 2 = E(^ f + V f )P+P n 

n 

H 3 = Y! M LP+Pn (7) 

n,m 

4 h^'M f (p+p+ + p p) 

1 1 2^ 1V1 nm \ r m r n ' r m*n) 

n,m 

The operators P n form the Paulionic algebra, i.e. they satisfy the following commu- 
tation relations. They anticommute like fermionic operators on the same site: 

PnPrt + PrtPn = 1 {P+f = = (8) 

and commute on the different sites like bosonic ones: 

P n P+ - P+P n = P+P+ - P+P+ = 0. (9) 
In the second quantization representation the operator of dipole momentum p n looks 

like 

Pn = VV (Kf b n0 + Ko b nf) 

or, in excitonic operators 

where is the dipole momentum of the 0-/ transition. In these terms the matrix 
elements can be rewritten as: 



Ip0/]2| 
Ifn 1 


n — m 


2 - 3(pV , (n - m)) 2 




n — m 


5 
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For the 1-D chain one can use the nearest neighbor approximation, so in the following we 
shall carry out summation in (|7]) over m — n ± 1 only: 



— E 

n,m n,m=n±l 

This is strongly different from the 3-D situation when one can not neglect distant mo- 
lecules, the number of which increases as third degree of the distance, and therefore the 
term Yl\ n -m\>i ■ ■ ■ therefore has the same order as the nearest neighbors contribution. 

One should note that the coefficient Mf nrn of the term P^Pn ls the same as that of 
the term P^P n (see (||, |) ) up to the effects of the intermolecular electron exchange 
which we will not take into account. This means, that the problem of the optical response 
is naturally resulted in XY-model rather than XX-model which was frequently treated 
earlier |,||. 

Now we will recall the formula for the linear response function. Let the aggregate be 
exposed in an external time dependent magnetic field. The energy of the aggregate field 
interaction is equal to that of the dipole momentum on the molecules in the electromag- 
netic field: 

H * = J2Pn E n{t) (P+ + P„) 
n 

(for every site the axis is chosen to be along with the molecular dipole momentum direc- 
tion). 

In linear on the external field amplitude approximation the most general relation 
between the statistical average of dipole momentum P n = (p°^(P+ + P n )) and the elec- 
tromagnetic field can be written as 



Pnit) = E / dT C: nm (T)E m (t - T) 

™ Jo 



where 



°n m {t) = Tr ([(P+ + P m ), (P+ + P n )(t))p) (10) 

is the linear response function, p = e _/3H / Tr e _/3H is the statistical operator, f3 is the 
inverse temperature. 

Making use the identity 

Tr (e* m (P+ + P m )e- lH '(P+ + P n ) P ) = Tr ((P+ + P m ) e -* m (P n + + P n )e iUt p) 

one can write down expression ([H]) in the form 

o- nm (t) = —T-, — {n nm {t) - n mn (-t)) (11) 
in 

where 



Ttr. 



l (t)=Tr((P+ + P n )(t)(P+ + P m )p 
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Expanding the last expression we get the final formula for the linear response function 
via the correlators of dipole momentum operators: 



Vnm(t) = ^ [(P+(t) P+> + <P+(t) P m ) + (P n (t) P+> + (P n (t) P m )- 
-{P+(-t) P+> - (P+(-t) P n ) - (P m (-t) P+) - (P m (-t) Pn) 

We can still simplify Eqn. (|l2l) . Indeed, due to the relation: 

Tr {e m P m e- im P n e'^) = Tr (e-^ H P+e iH *P+e" iH *) 

we have 



(12) 



(P m (t) Pn) = (P+ P+(t)) = (P+(-t) P+) 

Besides, the correlator (P^(t) Pn) is equal to (P m (r) Pn)\i T =a-it because of: 
Tr (e-^ n P+e- im P m ) = Tr (e- im P rn e~^ H P" 



Finally, using the time and the chain translation invariance we derive the following ex- 
pression for the response function a nm {t): 

°nm{t) = ^^(lm(P n + (t) P+) + Im(P+(-t) P+) + 

+Im(P n (i) P+) + Im(P n (r) P+}| . (13) 

So, we have got the expression for the linear response function cr nrn (t). Let us emphasize 
that there are two different kinds of correlators in (|i"3"|): (P„ f (t) P+) and (P n (t) Pm)- m 
the next section these correlators will be found for the general case of the cyclic XY-model. 



3 Correlators of XY-model 

In this section we will consider the calculation of correlators which are needed for the 
evaluation of the linear response following formula fll3|). 

The results of this section are valid for an arbitrary XY-Hamiltonian disregard to 
specific relations between its coefficients, so we can write again Hamiltonian ([3], ^) in the 
general notations: 

M M i M 

h = E^»+*E E ^» + ^E E (e+p+ + PmP») (w) 

n=l n=l m=n±l n=l m=n±l 

where the cyclic boundary conditions are imposed and P+, P n are paulionic creation 
and annihilation operators obeying the commutation relations (|8|, |^). For excitons the 
parameters are e = Ae/ + Vf, h = 7 = M£ m We have also eliminated the constant 
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summand H 1 ([|), because statistical characteristics are independent on a constant shift 
in the energy operator. 

It is convenient to start with the consideration of the partition function of the system. 
After that the correlators (Pi(t) P^ +1 (0)) and (Pi(t) P/ +1 (0)) will be evaluated. 

3.1 Partition function 

As usually, we make use the Jordan-Wigner transformation 

m — 1 m — 1 

, . E 4 c fe , , , , E 4 c fe 

P - _1 M C P+ = C +f — l)k=l 

to rewrite operator (0) in the fermionic Hamiltonian: 

M M M 

H = 6 C n C n + hJ2 ( C n C n+l + 4+l C n) + 7 ^ (44+1 + C„+lC n ) (15) 
n=l n=l n=l 

with the following boundary conditions for the fermionic creation and annihilation oper- 
ators: 

M 

CM+i = e-^ Cl c + +1 = c +e i7riV N = J24c n . (16) 

n=l 

Let us denote the orthogonal projection operators on states containing even and odd 
number of particles by P + and P_ respectively: 

P± = \ (l ± (-1)*) (17) 



and consider the operators H± of form (|15D with the boundary conditions 



cm+i = =Fci cjf +1 = =FCi . 
In the momentum representation 

c m = £ e ipm c p , A ± = { P e (-7T, tt] : e* Af = T l} 



M peA± 

the operators H± take the quasidiagonal form 

H ± = [e p CpC p + 7 sin p(c+c+ p + c_ p c p )| , e p = e + 2h cos(p) 

P eA ± 

Then the partition function Z with the chemical potential \i is represented as: 

Z = Tr {exp(-/3(H+ - /iiV))P + } + {Tr exp(-/5(H_ - //iV))P_} 
that can be rewritten in the form: 
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Z = - ^Tre _/9(H+ ^ ) + j te -fi( H --^) + 

+ Tr(-l) V^ H +-^ - Tr(-lfe-^—^) (18) 



Such form is convenient for the calculation of the partition function. Indeed, each term 
of Eqn. fll8|) is a partition function of a quadratic Hamiltonian which is quasidiagonal in 
the momentum representation. Moreover, the moments k and —k are only connected. 
So evaluating traces in each block and taking the product over all momenta we have the 
following result for the partition function H: 



1 
2 

T± ( ry± 



Z = -{Zf + ZJ + 1/Z+ - 1/Z b ) . (19) 



Here Zj (Z b ) are the fermionic (bosonic) partition functions 

Z± = -Q e (3E(k ± )/2 A + e -f3E(k ± y 

fc± 

Z± = e^*^ 2 JJ e -^( fc± )/ 2 (i - e - ^**)) _1 
fc± 

for systems with the energy spectra 

4 7 2 sin 2 (A; ± ) 



E(k ± ) = ( 



e k± — A* 



1 



\ (e*± - m) 2 
and the antiperiodic (periodic) boundary conditions: 

fc+ = ^(m + l/2), fc_ = ^m, m = 0,...,M-l. 

Using formula ( jl9|) one can calculate all equilibrium thermodynamical quantities. The 
partition function will be used further for calculations of the correlators. There we will 
put n = 0. 

3.2 Correlator (P±{t) P/ +1 (0)) 

In this subsection we concentrate our attention on the correlator (Pi(t) -Pl+i(0)). It will 
be calculated by the coherent states method adopted for the case of the cyclic XY-model 
in Ref.|P-l]| (the final result of the paper contains a mistake corrected in Ref.Q). 
By the definition, the correlator we are interested in is 

(Pi(t) = ^Tr (e-^ H Cl e- ffl *c+ +1 (-l)En=i^-«) (20) 

As in the previous section it allows us to rewrite correlator fl20|) in the following form: 
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(Pi(t) Pf+M) = ^ {K+ + K+ + Kj - K^} 



(21) 



where we have introduced the notations: 



Kf = Tr ^ (/3 ^ )H± Cl e^ H T*c+ +1 (-l)SLi c " c ^ 

: Tr (V (/3 -*' )H± Cl e- m ^ci +1 (- i)D^+i c - c ") 



Let us calculate accurately the first term Kf . We will consider the basis of coherent 
states 



p± 



with the following well-known properties: 



J P ± 



id 



-a 



p±i 



The variables £ are Grassmann numbers and are taken to be commuting with operators. 
In the Bargmann-Fock representation correlator Kf is expressed as follows 

Kf = I n <+ < / n < + < + / n < < >< 



•itH" 



S P+/ \ S P+ 



,0 \ 

Vi-/ 



(22) 



We will consider each multiplier in details and start with the third one. Making use the 
identity (N is the normal ordering symbol) 

e ac+c = Nexp ((e a - l)c+c) 
one can rewrite the operator (— l)S*=i C X C * in the normal ordered form: 



_l)ELi c fe c <= = Nexp< 



"2 E sVp' + <p' + 



p+.p+ 



with the matrix V* : 



1 sin 5 (p' + i.,, 



p+p '+ M sm§(p' + -;P+ 



e W+-P+)(L+l) _ 



(23) 
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Then, substituting the momentum representation for the operator c^ +1 

p-lTr/4 

4 +1 = ^Ee--+ (L+1) c+ (24) 
V JVi p+ 

we come to the expression for the third multiplier in fl22"|) : 



where 

<„v = - 2 < pV + W + • 

Let us now consider the first multiplier in the integral (|22|). Using ( p4"|) one can rewrite 
it as follows: 

f o +|e -(^ )H+Ci|e i + \ = ^=Ee ip "'^(g + |e-^)H + | _£\ (26) 



— 17T/4 



M v ,n 

To calculate the matrix element we have to present the operator e~^~*^ H+ in the normal 
ordered form. First of all we will extract terms in H + with the momentum p — tt (if 
it is presented at all). All other terms can be grouped in pairs with opposite momenta, 
— p + }. Then the exponent can be put down: 

e -Q3-tt)n+ = e -(/9-it)H» -q e -{p-it)n p + ^27) 
where the first multiplier is only presented if M is odd and 



B. p + = e p (c+Cp + c+ p c_ p ) + 27 ship (c_ p c p + c+c] 



One can derive the following expression for the normal ordered exponent of the last 
operator: 



e- (/3 ^ )H+ = e-(^ + K p+ (P - *t)Nexp j [ Kp+{ p_ it) - 1J (4 + c v + + ct p+ c. p+ ) + 
+a p+ (l3 - it) (c+ + c+ p+ + c_ p+Cp+ ) I 

where 

KM = cos 2| e ^ + sin 2 |e-^ 



1 



ct p(^) = 17Tff\ sin 6p^P E P 
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and 9 P is so-called Bogoliubov's angle: 



sin 6L 



27 sin p 
E n 



cos 6 P = , 



E = £ 

J-Jp Op 



1 



47 2 sin 2 p 



+ 1 . 



Now one can get the expression for the matrix element in (|26|): 



ll + e-(^ H + -e p+ ) = II e-V^K^ip - it) 



x exp< 



{p+>-p+}^i- 



k K p+ (P-it) ^ p+ ^ p+ 
-a p+ (f3 - it) + e p+ e p+ ) } exp (-e 



-(j3-it)eir FO tl 



This expression allows us to consider the product over all momenta p + (instead of the 
product over the pairs — £>+})• Then the terms with the momentum k = ir can be 
taken into account in the product over all p + and the term fl2"6|) can be presented as 



3 «7r/4 



M E eJpV 4 (ne- (Me ^ + (/? - it) 



x exp 



\ FO cl _ °P+ )P ^) ( cO cO 

' K p+ ({3 - ity p+ ^ p+ 2 



/to fo + *i a 



(28) 



In close analogue with this term we can derive an expression for the second matrix 
element in (|22|) . It is given by 



e^ H " Cl 



p+ 



Y[e- it£ *-K p _(it) 1 x 
,p- 



x exp 



v ^ p „(it)^-^- 2 

where the variables £ p _ are connected with the variables £ p+ by the unitary transformation: 



£p- ~~ 5Z^p-p+£p+ ' £p- ~~ 5Z^p+ U } 



p+p- 



p+ 



1 M 

P ~ P+ ~ Mj^ 



Substituting (|25|), (|28|) and (^) in the integral fl22|) one can get the final integral expres- 
sion: 
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2 



\p+ p- J M P+ ,p> 



X 
X 



/ d£° de° d^ 1 df d^ 2 exp {er - ^ - ? e 2 } x (30) 
ex P {-?D { ?e+ + e ^ - eB?e + i x u + D^ue + 

where we use the matrix notations: 



n(i) _ HiW i ^ Rf 1 ) _ A °"p^ ~ ffi 
D ~ dm S{ Kp{p _ it) ) , B vpl --6 P> - P , 



D{2) =diag(^y), = -5 P ,- P 



a p (it) 



indices "+" and "— " at matrices point the momentum space (A + or A_) in which 
these indices are considered. 

One can see that can be written in the same form as K^. There is the only 
difference in the third matrix element in ([22|). It leads to replacing of the matrix W(_ by 
the matrix W+. The matrix can be easy calculated and results in 

h 2 fl , . w . I sin|(p' -p+)(M-L) f 

W», = 5 P+P , - — exp -i(p' - P+ )(M + L + 1)\ 2 T + ! v } = -W* 



p + p' + »p + p + M 12 + ™ ' ' ! \ sini(p' + -p + ) " 

That is why we have the general expression for Kt and in (|30|). 

To obtain Kjr b s one should replace the matrix U by the matrix U + , "+" by "— " (and 
visa- versa) in Eqn.(pO|). It can be easy understood from the calculating procedure. Thus 
KJq>) * s P resen ted by: 



K m = (u^ u)£p - K P^-^U^ lt£p+K pA^)) ^ E e^'--^-( L+1 ) x 

\P— p+ / p-,p'_ 

x / de° df d^ 1 de 1 dp dp exp {e°e° - pe 1 - e 2 e 2 } x 

x exp {-^D W e+ + |° Si 1 ^ - ^fii 1 ^ 1 + t l UDfu + e + 

+t i uB®u t ? - euB®u + e - ?wi (b) ?} (31) 

So we have obtained all terms in ([H]). Now we should calculate the integrals in expressions 
([ft],[n]) . To do it we introduce four generating functions (fj, rj) and Ij^ (fj, rf) as follows: 
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it® (a, v) = J ae dp dp dp exp (er - pp - pp) x 

x exp {-^D^e + P4 1} p - pi^P + ?U + D™U£ 2 + 

+z i u + B™uz i - eu i B m ue - ?wi {b) e + m 1 + ?v} ■ 

^) = / d ^ d ^ d ? d ^ 2 ex p - ?e - ?e) x 

x exp {-?D^e + pfi a) p - p^p + p^i% + p+ (32) 
+?ub?v*? - euB { ?u+e - ?w m e + m 1 + e~M • 



The detaled calculation of these integrals and corresponding derivatives is placed in Ap- 

^± . 
L Ap- 



pendix A. Here we only put down the final results for Kf 



Kf m = (-^^^PfCGfVdet^x 



L /(6) 



> l n e~ (/3 ~ l ' )£p ± K p± ((3 - it) J] e~ ltep T X Pt (it) x 
x -L £ e^-^+D (gf ) + Ti (b) (g£ (b) ) "* 5 ± ) _1 (33) 



p±,p ± ■ p ± p± 



r 

with the following notations 



gf) = wWdP + u+dWr-U 

= 2{^ + /W i?i 1) (^ (b) ) t + f/ + fii 2) J R_c7} (34) 
S+ = -2 {s^ + U l B m R_u) 

Q f_(») = w m D W + UD^R + U + 

T m = 2! [ W m B W (w m ) t + UB^R + U^ (35) 
S- = -2{B { } ] + UB i + ) R + U + } 

and the matrix G±^ has the form 

G m_( s ± (Qf»)'\ 
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where the matrix R is defined as 

R = |(Z}( 2 )) 2 - 4 (B^) 2 

3.3 Correlator (P^{t) P/ +1 (0) 
Let us now consider another correlator 



diag ■ 



K 2 Mt) 



1 + K*(it)a 2 p (it) 



Pl+i(0)) = ^Tr fe- MH c+e^ tH c^i f-D^Li <£<* 



-L+1V 



In close analogue with (-Pi(t) one can rewrite the correlator as: 



where 



(37) 



(38) 



L + 



#± = Tr ^ e -(^-it)H ±c + e -itH Tc + +i (_ 1 ^ fc=1 



The calculating procedure is the same as in the previous subsection. But there are 
several differences on which we will concentrate here. 

First of all, let us consider expression Q22p. There is the only difference here for K'j 
in two first matrix elements in (j22|). For Kt they are: 



cfe~ ltH - 



t. 



p+ 



to provide the normal ordering. We have estimated all matrix element for Kt in the 
previous section in details. From these calculations one can see that for K~j we should 
change the matrix D by the matrix — D in (|30"D, change the sign in one of the exponents 
in front of the integral and put the multiplier e~ ln / 2 = —% due to the properties of Fourier 
transformation. 

So one can infer the expression for K^ b y. 



K 



f(b) 



J| e -(0-«)s, + Kp+ (p _ if ) Y[ Q- it£p ~ K p _ (it) 



H) 

M 



e -¥+-ip+(L+l) 



X 



r2 



< / d^° df d^ 1 d^ 1 d? de ex P - ?e - ?e} ^ 

< ex P {^°D^e+ + fs<_v - ^B^e - l l U + D {2) Ue + 



(39) 
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where all notations have been introduced in the previous subsection. Using the same 
arguments as above for Kj,^ we can derive the analogues expression for Kj, b s 



.p- p+ / p-,p_ 



x 



/ d^° d£° d^ 1 d^ 1 df d^ 2 ex P {er - ^ - e 2 x 
ex P {z°D w e+ + e^e - e 1 ^ 1 - ^uD^u+e + 

+l l UBfu t ? - £ 2 tLE?f t/+£ 2 - f iy/ (fe) e } (40) 

The next step is to calculate the integrals in (|39| , |40|) . We will consider just one of 
them and write another one by analogue. Let us introduce the generating function for 
the integral in K^ b y. 

4 {b \vi,v2) = J ae df d^ 1 d^ 1 df de ex P - ?e - ?e} * 
x ex P {^D^e + f - es^e - i i u + D m ue+ 
+i i u + B m ui i - eu i B m ue - ?wi {h) e + ev + e 2 ^} 

Integrating in the same way as in the previous section and neglecting terms with 771, 
771 and 772, rj2 one can obtain 

Jfafarh) = (-l)^ 1 Pf(Gf V det R ~- ex P {mFl {b) V 2 } 

Here matrices Q, S and T are the same as in fl5ip and matrices and are 

determined by 

d m = ( s ± -(g£>) 

= (2s£) + D^(wl {b) y (Qi^y 1 s- ± ) (g£ (fe) + T| (b) (g£ (fe) )~ ' s ± x 1 

Using this generating function we can write expressions for all K in (|38|) applying the 
same calculations as in the previous section. The result is given by the following formula: 



K f J b) = (-1)™ ne" (M)ep± ^ ± (^ - ^ ne" itep ^ PT («t) 

\P± PT , 

x P^G^Vdeti?" 1 ^ £ e -y ± -, P± (L + i) (Fi (fc) ) p±p/ (41) 

p 1 y. 
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Here and matrices >S_ and T_ are the same as in (35). So formula (41) gives us all 
the necessary terms to calculate correlator (|38|). 



4 Correlation functions for 7 = (XX-model) 

Up to day there are three original works devoted to the calculation of the correlators for 
the XX-models at finite temperature. The first of them is the paper by Colomo, Izergin, 
Korepin and Tognetti |7| where the correlator for the XX-model was calculated in the 
thermodynamic limit. The correlators for finite cyclic chains were calculated by Suzuura, 
Tokihiro and Ohta in |J and by Ilinski, Kalinin and Kapitonov H for the qXX-model 
(which produces the XX-model for the particular value of the deformation parameter 
q = —1). Thus, there are three crucial points for the verification of the results obtained 
in this paper. Besides this, it is possible to compare our result with that of calculated by 
the definition for small chains. It means to represent the Hamiltonian and the creation- 
annihilation operators as matrices 2 M x 2 M for the M-site chain and then to evaluate the 
trace of the matrix e~^~ lt ^ H Pie~ ltn P^ +1 . The result of the comparison will be reviewed 
in this section. 

Let us start with the finite XX-chain correlator. One can extract it from the formulae 
(0-0) putting 7 = 0. Then 9 P = for all p and E p = e p . 

{Pi® PUM) = 7^(Kl + KL + K b + - K b _) (42) 



K Kh) = ^_-QAf(Af-l)/2 JJ e -itep pf f (Fl ( b) Y\ ^(pf®- 1 

P eA+ \ ~P± / 



G) 



P eA ± 

= e~ 4 ' (£)M detFl {b) Tr(Fl (b) G) (43) 
where 



,/(&) 



,/(*>) 



w m D {i) 



U+D^R.U 



UD^R+W 



7=0 
7=0 



G 



p.p 



M 



_ e -ip(L+l)+ip' 



(44) 



If we rewrite this formula in the coordinate representation we will get the formula for the 
correlator in the qXX- model obtained in [[| for the case q = — 1. 



^det mm$+i 



(45) 
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( pf\ , - J_ V P iter+ir(m-m') , J_ pip{m-m')+(it-0) £p ( \ _ r, J 1 if TO < 

^± w - M 2^ e + M e l v i ^ x ioifm> by 

(f?b\ , J_ pite r +ir(m-m') L "S^ p ip(m- m ')+(it-l3)e p f i _o v I 1 if m < 

One should note that our result is different from that announced in || . First of all in 
our opinion there is a misprint in signs of two last terms in Eqn.(15) of Ref. |J. Moreover, 
the numerical calculations for t — and particular values of parameters ({e — l,h — 
5}, {e = 1, h — 2}) reveal the difference between the results of and the calculation by 
the definition for both choices of the sign. 



Eqs.( f42| , f45|) coincide with the correlator calculated by the definition for the number 
of sites M = 3. In this case it is possible to compare the analytical expressions for the 
correlators (the calculation was performed using the system of symbolical calculations 
REDUCE). 

Let us now consider the thermodynamic limit (M — ► oo) for the XX- model (7 = 0). 
At first we will treat the partition function. The "bosonic" contributions are equal 
and, hence, disappear. By the same reason the "fermionic" terms are equal and cancel 
the factor 2 in the formula (p!9|). This is due to the absence of the difference of momentum 
spectrum A + and A_ (in the thermodynamic limit the boundary conditions are irrelevant). 
So the partition function of the XX-model in the thermodynamic limit is equal to the 
partition function of free fermions on the infinite chain. 

In the same way the correlator for the XX-model in the thermodynamic limit is ob- 
tained from Eqn.(fi^) putting K+ = KL = K, = K h _. So in the thermodynamic limit 
we have 

(Po(t) P L + (0)> = § (48) 

where 

Z = + e-^f) (49) 

v 

is the fermionic partition function and 

K = e~^ M det(F)-L £ e *L F -i , (so) 

F = R(it -(3) + R(it) - 2VR(it - (3) , 

R(t) = diag{e T£p } 

and the matrix V is given by Eqn.fl2"3p. 

It is easy to show that the matrix F can be represented in the following form: 

F = 0- 1/2 [l - 26 1/2 V6 1/2 ]6- 1/2 R(it - (3) 

where 9 is the following diagonal matrix 



9 = diag (- — ^—5—1 
1 1 + e? e p J 
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So we obtain 

det(F) = Ze it£M det[l - 26 1/2 V6 1/2 } 

hence 

(P (t) P£ (0)> = det[l - x 



X M , 

Let us consider the following unitary transformation of V: 

V = U~ x vU 

where 

rj = dia g {e^ L+1 )/ 2 } « ■ - 2 



sin^ 



Then the correlator takes the form 



(P (i) (0)) = det[l - 26 1 / 2 v6 1 / 2 }-^ £ e i(p+p ' )L/2+i(p - p ' )/2 iV P '(/3 - it) x 



x^4 2 ,[i-2^w/ 2 ] p -% 2 

det[l - 2^1/2] _ 25,1^-1^ 



_ . x 

p.p 



xe *(p+P / )V2i[ e »(p-P / )/2+(^-it)v _|_ e i(p'-p)/2+(/3-it) £p j 

And finally we obtain the following expression for the correlator 

(P (t) PjJ-(O)) = det[l - 2# 1/2 i# /2 ] Tr {# 1/2 r(/3 - zt)6» 1/2 [l - 26 l/ \6 1/2 \- 1 ) 
d 



, det[l - 20 1/2 v9 1/2 + z9 1/2 r{(3 - it)9 1/2 } 
oz 



2 = 



where v is defined in (1511) and 



rpp ,( r ) = e *(P+P / )i/2i[ e i(p-p')/2+r£ p , + e i(p'-p)/2+r ep j 

In the same way expression for correlator (P + (t) Pl{0)) is obtained 

(P+(t) P L (0)} = det[l - 2# 1/2 ^ 1/2 ] Tr {# 1/2 r(^)# 1/2 [l - 26 1/2 v6 1/2 }- 1 } 
^ det[l - 20 1 / 2 ^ 1 / 2 + zd^rme 1 ' 2 } 

9z z=0 
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The form of expression (54) is similar to that in Ref. 0. Let us substitute t = in 
Eqn.fl): 



(P+ p L ) = det[l-26 l / 2 v6 1 / 2 }Tv{6 1 / 2 r6 l / 2 [l-26 1 / 2 v6 1 / 2 ]- 1 } 
d 



dz 



det\l-29 1 / 2 v9 1 / 2 + z9 1 / 2 r6 1 / 2 } 



2=0 



where 



r PP > = tw(0) = e 



i{p+p')L/2 



cos 



(P - P') 



(55) 



(56) 



Now we can see that expression Q55D is different from Eqn.(39) of Ref. [[?]] in two aspects. 
At first, the matrix r contains the additional factor cos . Secondary, in the matrix u 
we should replace L by L — 1. We will show that these two expressions for the correlator 
(P + Pl) are the same. Indeed, we can write down from the definition 



W Pl) 



Tr [e-^Po+P^j 

d£o d«^o d^i df i exp[£ £o - x 

x(eo|e^ H iei)(eiiNc+ex P {-2x: c + Cfc } Ci ieo) 



(57) 



k=l 



Using the method described in Section 3 we obtain the following result 

(P+ P L ) = det[i-2e 1 / 2 ve 1 / 2 }Tr{e 1 / 2 fe^ 2 [i-2e 1 / 2 ve 1 / 2 ]- 1 } 

d 



dz 



where 



det[l- 2e x W + ze l ' 2 r9 l l 2 } 



f , _ p i(p+p')L/2 
PP 



(58) 



2=0 



Vpp' 



sin 2 (p' — p)(L — 1) 
sin | (p' — p) 



As CqCl = Cq exp[iTrcQC }cL we can rewrite ([57]) in the following form 
(P + Pl) = \\ d£o dfo dfi dCi exp[£ £ - £i£i] * 



x(eo|e^ H ie 1 )(eiiNc+ex P {-2x; 1 c + Cfc } Ci |eo) 

k=0 



(59) 



Repeating the above calculation we can come to expression fl55|). It follows from the 
previous consideration that expressions fl5"5"D and (^) are in fact the same. 

So it is not now surprising that expression (154]) for correlator (P + (t) P^(0)) is different 
from that in Ref. [7J because they are derived with different methods and as we have seen 
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above it is difficult to transform one form to another. However, due to the verification 
for finite M and for t = 0, M — > oo we can be sure in the correctness of our formulae in 
M — > oo limit and t ^ 0. Finally, our result is obviously simpler. Indeed, expressions 
(0,0) for the correlators contain the time dependence only in the matrix r and trivially 
generalize the t = case. 



5 Conclusion 

Now we can summarize the results of the previous sections. Formula (|13|) together with 
formulae (21], |33], |38|, [41]) for the correlation functions give the final expression for the 
linear response of the molecular aggregate. In our opinion, the result has the following 
advantages: 

1. It can be used for the numerical calculations to investigate the finite size effects 
and the closure effects for cyclic aggregates. As we could note during computer 
calculations for the XX-model, to deal with our expressions is significantly easier 
than with the corresponding formulae of Ref . |§ 

2. We can say the same about the resulting representations for the correlators in the 
thermodynamic limit. We cannot check our expressions directly in this case but the 
positive results for small finite chains assures us. 

Besides, we want to add that our Fredholm determinant representation of two-point corre- 
lation functions of the XY-model also allows to write differential equations and to calculate 



their asymptotic by the standard way [12 



In conclusion, in the paper we considered the linear response function of a gas of 
cyclic aggregates in the framework of Frenkel excitons. To do this we reminded that 
the problem is effectively equivalent to the calculation of the correlation functions for 
the cyclic 1-D XY-model. We obtained the exact formulae for two-point correlation 
functions at finite temperature for the model and compared them with previous results. 
In particular, we showed that for the XX-model our result for the correlation function 
< P n (t)P+ > coincides with the straightforward calculation (using the definition of the 
correlation function) for the simplest nontrivial case (M = 3) and is reduced in general 
to the result of paper || for an arbitrary number of sites. We also found that the result 
of Ref.|| for the correlation function is incorrect. The thermodynamical limit of our 
formulae was discussed. 
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00548 and partially (K.I.) by the UK EPSRC under Grant number GR/J35221, Euler 
stipend of German Mathematical Society and grant INTAS-939, (V.M.) by President 
stipend, (A.S.) ICFPM Fellowship, ISF grant No A602-F. 
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A Appendix 

In this appendix we will calculate integrals ( |52"D using the following formula: 

J d(f> exp (~<f>K<f> + <f>aj = Pf(K) exp (-aK^aj (60) 
where </> and a are Grassmann vectors and Pf(K) is a Pfaffian of matrix K. 



Let us consider 1^ and integrate over £ , £ first, using formula (|6"0~D , where we will 

put 

*-($) 

Then the result of the integration is given by: 

exp {fWf 6) 4 1} K^)*? - FwIVd?? 

Now we will integrate the remaining expression over £ 2 , £ x to keep aside the terms with 
77, 77. The result takes the form: 

'deti^exp {?U + B®IL.U? - £ 1 U t B®IL.UZ 1 - fCZ+p^)- 1 ^^ 1 } 



where for a convenience we have introduced the diagonal matrix R (|37). The third 
integration over £ 1 , £ 2 gives us the final result. We will neglect the terms in the exponent 
containing fj, fj or 77, 77 because they will vanish if we calculate the integrals in Klfn* using 
generating function 7^(77,77) (by taking derivatives on r\ and 77). So we obtain the 
following expression: 



iUv,v) = (-i)^Pf(Gf>) V / d^exp|77 (gf > (gf>)-*s + )~ , 

where 

G m = ( s + («?>)' 

and for the sake of brevity we have introduced the new matrices Q+ b \ T+ and S+ which 



are determined by (|34|). 

In the same manner we can present a result for the generating function Ij^ : 



i m (v,v) = (-i) M ^ ii Pf(^)v / d^?ex P 77 (ql® +t ™ (Q m y ts -T V 
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Here is defined by changing all indices "+" at operators by indices "— " in G+ and 

Q_ ,T^ b \ S- are determined in (|35|). 

To obtain an expression for the integrals in Kj 1 ,^ one should consider 



d%' + dr] p+ 

and then the final result for KjT/ b \ can be written as: 



J] e- {f3 - lt)£ ^ K p± {p - it) n e~ ite ^ K PT {it) 

yP± PT 

J_ v e ^ ± (w)U) + T /W(^))^r 



p±,p' ± x ' p ± p± 
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